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We performed extensive Monte Carlo simulations of the irreversible adsorption of polydispersed 
disks inside the cells of a patterned substrate. The model captures relevant features of the irreversible 
adsorption of spherical colloidal particles on patterned substrates. The pattern consists of (equal) 
square cells, where adsorption can take place, centered at the vertices of a square lattice. Two 
independent, dimensionless parameters are required to control the geometry of the pattern, namely, 
the cell size and cell-cell distance, measured in terms of the average particle diameter. However, to 
describe the phase diagram, two additional dimensionless parameters, the minimum and maximum 
particle radii are also required. We find that the transition between any two adjacent regions 
of the phase diagram solely depends on the largest and smallest particle sizes, but not on the 
shape of the distribution function of the radii. We consider size dispersions up-to 20% of the 
average radius using a physically motivated truncated Gaussian-size distribution, and focus on the 
regime where adsorbing particles do not interact with those previously adsorbed on neighboring 
cells to characterize the jammed state structure. The study generalizes previous exact relations 
on monodisperse particles to account for size dispersion. Due to the presence of the pattern, the 
coverage shows a non-monotonic dependence on the cell size. The pattern also affects the radius of 
adsorbed particles, where one observes preferential adsorption of smaller radii particularly at high 
polydispersity. 

PACS numbers: 02.50.-r, 68.43.Mn, 05.10.Ln, 05.70.Ln 



I. INTRODUCTION 

Monolayer colloidal films grown by particle adsorption 
on substrates have a long standing interest from both 
scientific and technological perspectives. Recently, the 
interest has shifted towards structured films at the sub- 
micron scale with the intention to achieve the nanoscale. 
From the experimental perspective, the emphasis has 
been on photonic crystals, quantum dots, heterogeneous 
catalysts, and microarrays [IHZ]. In particular, there has 
been progress on achieving highly reproducible and reg- 
ular submicron patterns jil, 5|i '81412). and control of the 
final structure of such films is of paramount importance 
[T31H4]. In this context, new probing techniques to exper- 
imentally follow the kinetics of the films along with novel 
theoretical methodologies are required. Specifically, the- 
oretical studies on the influence of a pattern on the ir- 
reversible adsorption of particles have so far remained 
restricted to the monodisperse case [TMT?] . Studies on 
polydispersed particles have been performed on regular, 
i.e., non-patterned substrates [18H22) . Nonetheless, ex- 
perimentally produced colloidal particles are inherently 
polydisperse, so there is a need to understand the effect 
of size dispersion on the morphology of the film. 
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We partially fill the gap by presenting a study con- 
cerned with the geometrical parameters and properties 
of the deposition process, and specifically on the inter- 
play of the adsorption of size dispersed particles on pat- 
terned substrates. The pattern consists of an array of 
cells centered on the vertices of a square lattice. Parti- 
cles are allowed to deposit inside the cells, which we con- 
sider in the present study to be squares of equal size. To 
this end, particles are considered to attempt adsorption 
on the substrate with the same probability regardless of 
their position, so that their motion is not, for example, 
affected by the hydrodynamics of the solution [231429] . 
Particles, which we model as disks, successfully attach to 
the substrate when their geometrical centers fall within 
a cell and also do not overlap previously adsorbed ones. 
If a particle overlaps a previously adsorbed particle or 
if its geometrical center does not fall within the cell re- 
gion, it is assumed that the deposition attempt fails and 
the particle moves away from the substrate. Moreover, 
once adsorbed they cannot either diffuse on or detach 
from the substrate [3D], so we consider the ideal case of 
irreversible adsorption. In this context of accounting for 
the basic (geometrical) parameters, the adsorption pro- 
cess is well described by the random sequential adsorp- 
tion (RSA) model [TH [T71 l5TH5o] . Landing particles are 
considered to have a truncated Gaussian-size distribution 
with dispersions up to 20% of the average radius as this 
represents a closer description of experimental systems. 
Finally, we focus our study on the regime where the ki- 
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netics within a cell is decoupled from the neighboring 
ones, i.e., the cell-cell distance is equal or larger than the 
diameter of the largest particle. 

The paper is organized as follows: in the next section 
we present the motivation for the present study along 
with a description of the model. In Sec. |III| we generalize 
concepts previously defined for monodisperse to polydis- 
perse particles, while in Sec. |IV| we present Monte Carlo 
simulations characterizing the film morphology. Finally, 
in Sec. IVlwe present our concluding remarks. 



II. MOTIVATION 

We address the study of the irreversible adsorption 
of polydispcrsc colloidal particles on patterned surfaces. 
The pattern takes a particular shape repeated over the 
surface and in the present study we use as a template 
the vertices of a square lattice. We further assume that 
adsorption of particles takes place solely within these re- 
gions we term as cells. We adopt the case of square 
cells for simplicity, but depending on the experimental 
technique used to generate these, they can assume differ- 
ent shapes with stripes and circles being experimentally 
produced [T^J 133 EE]. Though we define the pattern 
as regular and deterministic, one might consider relax- 
ing such conditions and consider random patterns. For 
example, an experimental realization of such a pattern 
can consist of three-dimensional hexagonal inverted pyra- 
mids on In-containing and Ill-nitride substrates, where 
colloidal nanoparticles of linear dimensions in the range 
5 — 30nm pack inside them |39) . Another example is the 
independent adsorption of colloidal particles of different 

sizes mam. 

The classical linear dimensions of colloidal particles are 
few microns down to about a micron, but recently the 
trend has been to use particles of sizes well into the sub- 
micron length scale. Similarly, the typical geometrical 
cell length has fallen below the submicron, and again, 
with a goal to approach the nanoscale. Even taking into 
account such experimental advances, the size of the col- 
loidal particles remains much larger than the length scale 
between binding sites, which are typically of the order of 
the linear dimensions of the substrate molecules. In this 
context, adsorption at the length scales of the cell linear 
dimensions can be regarded as an intrinsically off-lattice 
process. We assume that diffusion and detachment do 
not occur on a time scale commensurate with experimen- 
tal observations, so it seems reasonable to adopt the ideal 
case of irreversible adsorption [42 - 44] . We are primarily 
concerned with the effects provided by basic geometrical 
features of the pattern, and consequently discard other 
(possibly relevant) interactions beyond the excluded vol- 
ume interaction. Since particles only adsorb inside cells 
of the substrate, and do not interact between themselves 
(beyond the excluded volume interaction), a monolayer 
film is obtained [T^l l5TH5rjl |4"5H4"7] and the model is char- 
acterized by an asymptotic jammed limit where no more 
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FIG. 1. Generalized version of the phase diagram for cell 
linear dimension a vs. inter-cell distance /3. Cell-cell separa- 
tion P < /3 C yields the interacting cell-cell adsorption (ICCA) 
regime, while P > f3 c > 1 gives the noninteracting cell-cell 
adsorption (NICCA) regime. For cell sizes a < ct c < \f2/2 
we find the single-particle-per-cell adsorption (SPCA) regime, 
while for a > ct c we find the multiparticle-per-cell adsorption 
(MPCA) regime. 



adsorption events can possibly take place, and conse- 
quently a limiting value of the coverage is attained, i.e., 
the jamming coverage. 

There is presently substantial literature on the adsorp- 
tion of monodisperse colloids on regular substrates (i.e., 
without the presence of a pattern) [5TH551 1551 I5B1 28]. 
Here, we are primarily interested in addressing the ef- 
fect of size dispersion on patterned substrates. Unfortu- 
nately, even such an extension poses several challenges 
that we partially address. As described above, the pat- 
terned substrate is fiat with particles being allowed to 
adsorb inside the cell regions with a uniform random dis- 
tribution. This is clearly an idealization as one expects 
different adsorption characteristics at the cell boundaries, 
but we do not address such effect in the present study. 
Furthermore, we take the colloidal particles to be disks. 
Strictly speaking, the particles cannot be modeled with 
full accuracy as disks since larger spherical particles can 
accommodate smaller spheres underneath them, and this 
effect was reported in the literature [49]. Besides this ef- 
fect, spheres of different radii have different overlapping 
rules than disks and even the meaning of what one con- 
siders as a cell, i.e., those regions of the substrate where 
adsorption can take place, becomes ambiguous. As an 
example of the latter, a pattern of cells in high relief 
would have different effective adsorption areas than one 
sunken relative to the interstitial cell space. In the case 
of sunken cells, particles of different radii attempting ad- 
sorption have different effective areas for adsorption. 



III. SOME EXACT RELATIONS 

We start by properly characterizing the various regions 
of a phase diagram as shown in Fig. [I] The effect of a 
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FIG. 2. (Color online) Particles attempting adsorption have a truncated Gaussian-size distribution, as described in the text. 
The two snapshots represent 20 x 20 unit cells region from larger simulations showing a typical configuration at the jammed 
state with a = 0.96. Cell-cell separation is large enough for the simulations to take place in the NICCA regime: At a size 
dispersion of l%, part (a), /3 = 1.5 > fi c — 1.02, while for a size dispersion of 20%, part (b), ft — 1.5 > /3 e = 1-4. In (a), 
only monomers, dimers, and trimers within cells are possible (in fact, trimers are not observed in this snapshot), while in (b) 
tetramers and pentamers are also observed. 



pattern has been previously studied in the context of the 
adsorption of monodispersed particles |15l I17j . As in the 
monodispersed particle case, we define two dimensionless 
parameters, namely, a, as the linear dimension of the cell, 
and, /3, as the distance between adjacent cells, both in 
units of the mean particle diameter, 2(r). In terms of 
the inter-cell distance j3, the diagram is divided in two 
distinct regions, namely, interacting cell-cell adsorption 
(ICCA) and noninteracting cell-cell adsorption (NICCA). 
In the former, adsorption of a particle at a given cell can 
be prevented due to overlap with a previously adsorbed 
one at a neighboring cell, while in the latter case, the 
distance between cells is too large for such an overlap 
condition to happen. Regarding the cell size, the value 
of a also splits the diagram in two regions, namely, one 
region where a single particle can be adsorbed, which 
we term single-particle-per-cell adsorption (SPCA), from 
another, where more than one particle can be adsorbed, 
which we term multiparticle-per-cell adsorption (MPCA). 
We are primarily interested in a range of cell sizes where 
the number of particles adsorbed per cell is one or, at 
most, a small number (less than six), since this provides 
the highest control on the actual number of particles per 
cell and on their sizes. 



A transition from ICCA to NICCA occurs when the 
distance between cells, /3 C , is large enough to prevent the 
largest particles, i.e., those of radius r max , from overlap- 



ping in the region between cells, i.e., 



(1) 



Notice that /? c > 1 makes the transition to NICCA at a 
larger inter-cell distance than that of the monodisperse 
case. Now, the transition from SPCA to MPCA occurs 
when the linear dimension of the cell allows a second 
particle to adsorb in a close-packed configuration, which 
depends on the radius of the smallest particle, r min . The 
critical value is 



(2) 



with a c < v2/2, i.e., at a cell size smaller than the 
monodisperse case. Only in the monodisperse limit one 
has r min = r max = r = (r) and previously reported critical 
values are recovered [T7] . 

Notice that the values of j3 c and a c separating the var- 
ious regions of the phase diagram are independent of the 
actual functional dependence of the distribution function 
in terms of the particle radii. In fact, they only depend 
on the extreme values of the radii (minimum and maxi- 
mum) as shown in Eqs. and ^ and the value of size 
dispersion. Now, if the limiting values of the radii are 

m 



(r) - 2a 



(3) 



4 



and 



(r) + 2cr, 



(4) 



then from Eqs. Q-Q, the critical values a c and (3 C can 
be put in terms of the dimensionless size dispersion by 
defining the latter in terms of the mean radius as 7 = 
er/(r), so that 



and 



a. = ^(1-27) 



Pc = 1 + 2 7 , 



(5) 



(6) 



which allows a useful re-interpretation of Eqs. ([l} and ([2]), 
apart from also reducing the number of independent pa- 
rameters to one. From Eqs. ^ and Q one has 7 < 1/2 
|51j . again due to the functional dependence adopted for 
r min and r max (Eqs. ^ and Q, respectively) in terms of 
a. Now, taking two distinct distribution functions of the 
radii, but with the same size dispersion, Eqs. ^ and ^ 
show that their phase diagrams are equivalent. 

In the present work, to systematize the effect of size 
dispersion on the adsorption process, attention is fo- 
cused on the regime where cell-cell separation, /3, is large 
enough so that adsorption on a given cell is not af- 
fected by particles previously adsorbed on neighboring 
cells since j3 > j3 c (NICCA regime). Depending on the 
value of a and the level of polydispersity of the particles, 
one can have one or more particles inside each cell. The 
effect of particle polydispersity is illustrated in Fig. [2] by 
snapshots of the jammed state for a = 0.96 for two val- 
ues of the size dispersion, namely, 7=1% (Fig. [2](a)) 
and 7 = 20% (Fig.gb)) [52]. For the remainder of the 
text we refer to values of the size dispersion in terms of 
their value relative to the mean radius. We denote as 
aggregates the set of adsorbed particles in a cell, and 
also name aggregates with a specific number of adsorbed 
particles, say, 1, 2, 3, . . . as monomers, dimers, trimers, 
. . . As in the monodisperse case, the inter-cell kinetics 
decouples so that one can follow the kinetics within each 
cell [17] . On the basis of this property we were able to ex- 
tend exact relations obtained for the monodisperse case 
to account for size dispersion in both the single and mul- 
tiparticle regimes. For example, in the SPCA regime, 
i.e., for a < a c , the mean adsorbed radius becomes inde- 
pendent of both the cell size and the actual distribution 
function of the radii, which equals that of the particles 
attempting adsorption. Consequently, the distribution of 
the radii of adsorbed particles follows a truncated Gaus- 
sian with the value of the dispersion also following that of 
the particles attempting adsorption. This can be better 
understood by considering the density distribution func- 
tion of the normalized radii in the film, pl^dr], where 
the normalized radius rj is defined as 




FIG. 3. (Color online) Adsorbed radii distribution for pat- 
terned surface with a equal to 0.4 (solid, red line) and 1.7 
(dashed, green line), and the bulk case (a — > 00) (dotted, 
blue line). For size dispersion: a) 1%, b) 5%, c) 10%, d) 15%, 
and e) 20%. The vertical lines represent the mean normalized 
radius, (77), of adsorbed particles. Notice that for a = 0.4 the 
corresponding value is 1/2, regardless of the size dispersion. 



As defined, p(rj) represents the fraction of disks of nor- 
malized radius in the interval ]rj, rj + dr/[ such that 



p(r/)dr/ = 1. 



(8) 



n 



(7) 



In the SPCA regime, since the radii distribution of ad- 
sorbed particles follows that of the particles attempting 
adsorption, the mean normalized radius of the truncated 
Gaussian is {rj) = 1/2. Figure [3] illustrates these argu- 
ments in the particular case of a = 0.4. In fact, it is 
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possible to obtain an exact relation for the coverage in 
the SPCA regime, i.e., for a < a c , as 



7T (1 



■7?) 



(9) 



where 7 T = a T /(r), with er T the effective size dispersion 
of particles attempting adsorption, i.e., <j\ — (r 2 ) — (r) 2 . 
In deriving the above expression we considered the area 
occupied by a set of polydispersed particles, with only 
one particle per cell, relative to the area occupied by a 
monodisperse set of particles. Equation ^ shows that 
in the polydisperse case the coverage is higher and that 
it increases by a factor of 1 + 7 2 relative to the monodis- 
perse case. Moreover, the coverage in the SPCA regime 
does not depend on the shape of the distribution func- 
tion of the particle sizes. As an illustration, in Fig. |4^a) 
the simulation results, to be detailed in the next section, 
show the initial decay of the coverage proportional to 
(a+/3)~ 2 ; deviations from this behavior coincide with the 
appearance of dimers (and larger aggregates), as shown 
in Fig. |4|b) . The solid black line in Fig. 4|a) reproduces 
the monodisperse case found in Ref. [17 (vertically dis- 
placed for clarity). 

Though it is not, in general, possible to calculate the 
jamming coverage for any value of the parameters a, /?, 
and 7 in the MPCA and NICCA regime, it is possible to 
relate the coverage to different values of j3 > j3 c by using 



gj(a,/8,7) _ f u + Po 



(10) 



where /3 a is the simulated value. The uncertainty for an 
ensemble of N samples of the coverage is given by 



\ 



N 



'ji 



N 



and using Eq. ( 10 ) one finds 



a 8 j{a,/3 ,j) \a + /3 



(11) 



(12) 



Equation ( 12 ) shows that the uncertainty of a particular 
simulation can be adjusted to the new j3 value. Equa- 
tions (|9])-(12) extend similar equations for monodisperse 
particles |17j to the polydisperse particle case. In the fol- 
lowing section, we analyze further results from extensive 
Monte Carlo simulations. 
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FIG. 4. (Color online) For the noninteracting cell-cell ad- 
sorption regime, at the jammed state, as a function of the cell 
size: a) Coverage (for f3 — 1.5); b) Mean number of particles 
per cell; c) Mean adsorbed radius. Each line type corresponds 
to different values of the size dispersion, namely, solid (red) to 
1%, dash (green) to 5%, dot (blue) to 10%, fine dot (purple) to 
15%, and dash-dot (light blue) to 20%. In (a) the solid (black) 
line reproduces the monodisperse case in Ref. [17] (vertically 
displaced for clarity). In the inset of (a) is magnified the re- 
gion 0.8 < a < 1.2 for 5% (lower line) and 10% (upper line) 
of size dispersion (details in the text). The error bars are 
smaller than the line thickness. 



IV. MONTE CARLO SIMULATIONS 

We provide some motivation for a truncated exponen- 
tial distribution of the size of particles attempting ad- 
sorption on a substrate, for which we foresee two possible 
scenarios. On one hand one can consider that the typ- 
ical size dispersions observed in experiments are Gaus- 
sian distributed, and the fact that one does not observe 



all the possible sizes is mainly due to processing and the 
timescale taken for the observation. From a theoreti- 
cal perspective, such approach poses difficulties on how 
to define the limiting coverage. Strictly speaking, for 
a Gaussian distribution one expects the coverage to ap- 
proach unity, since particles of any size are possible, so 
the system does not jam, but this is not physically ex- 
pected. On the other hand, taking a truncated exponen- 
tial as the distribution function of the particle sizes has 



G 



several benefits: it naturally accounts for the absence of 
the extreme values of the linear dimensions of the par- 
ticles effectively present [HI [20]) and removes the non- 
physical result of reaching a fully covered substrate. We 
define the truncated Gaussian-size distribution, P(r)dr, 
as 



P(r) = 



Aexp 




r min < r < r n 
otherwise 



(13) 



where (r) is the mean radius, a is the dispersion as mea- 
sured from a Gaussian distribution, and A the normal- 
ization factor. All particles with radius below r min and 
above r max are discarded on the basis of being extreme 
values not actually observed within a typical experimen- 
tal timescale [53] . The dispersions of the truncated Gaus- 
sian, cr T Gi and of the Gaussian, er, are related by 



<T TG 

a 



2 3/2 



0Fe 2 erf(V2) 



1/2 



(14) 



However, the actual cutoff values for r min and r max are 
given by Eqs. (|3| and Q, using the Gaussian distribution 
value, instead of the truncated Gaussian value, for a more 
direct interpretation. Consequently, the coverage in the 
SPCA regime is given by 



7T(1+/V) 

4(a + /3) 2 ' 



(15) 



To study the influence of cell size on various quanti- 
ties such as the coverage, mean adsorbed radius, den- 
sity of monomers, dimers, etc., and mean number of ad- 
sorbed particles per cell, we resort to extensive Monte 
Carlo simulations. From the experimental view, size dis- 
persions above 20% are rarely observed or of interest, so 
we consider size dispersions up to this value in the set 
of simulations [20]. As mentioned in Sec. [TTJ the tran- 
sition from single to multiparticle adsorption occurs at 
a c = (1 - 27) v / 2/2 (Eq. We also mention that 

/? = 1.5 > /3 C = 1.4 for the case of a size dispersion of 
20%. Each simulation was carried out on a substrate of 
500 x 500 unit cells and for an ensemble of 10 2 samples. 

We used an efficient algorithm, to be described in 
greater detail in a separate publication [54], while in the 
present work we present only an outline of its main as- 
pects. We note that similar algorithms exist in the litera- 
ture [T5 1 [T7 | H5 1 l2n i l5lH55]. The substrate is divided into 
a homogeneous mesh of squares and their size set to the 
minimum particle radius, r min . We denote the squares 
as mesh cells to avoid confusion with the physically rel- 
evant cells of the pattern. Each mesh cell is classified as 
empty, occupied, or shadowed and only empty mesh cells 
are tested for adsorption. Occupied mesh cells contain 
the geometrical center of a particle and, therefore, they 
cannot be available for adsorption. Similarly, shadowed 
mesh-cells are those where adsorption cannot take place 
due to the excluded volume interaction of a previously 
adsorbed particle. Once the number of empty mesh-cells 



falls bellow a critical percentage, typically below 2% of 
its initial value, the linear dimensions of the mesh cell 
are halved. Since mesh cells classified as empty can have 
part of their area shadowed, the fraction of shadowed 
area present in an empty mesh cell can be reduced by 
halving its linear dimensions and eliminating those off- 
springs that are now shadowed. On average, the pro- 
cedure reduces the total shadowed area present in the 
empty mesh cells. Such a reduction of shadowed area 
leads to a higher probability of acceptance of a particle 
attempting adsorption, and, consequently, to improved 
algorithmic efficiency as compared to prior of the halving 
procedure. Further details of the algorithm are discussed 
in Ref. [54]. 

We start the analysis by considering the density dis- 
tribution function of the normalized radii in the film, 
p{vi)drj. At values of a < a c the distribution of adsorbed 
radius follows that of the particles attempting adsorp- 
tion. It must, thus, correspond to the truncated Gaussian 
distribution, and the observed distributions for a = 0.4 
in Fig. [3ja)-(e) do corroborate the expectation. This is 
no longer the case for a = 1.7 where several particles 
can fit inside a cell, but the smallest particle sizes, on 
average, have higher probability of adsorbing than larger 
ones. Though valid at all values of the size dispersion, 
this effect becomes more striking for the highest values 
of the size dispersion. As the values of the size disper- 
sion increase p(r]) becomes more skewed towards smaller 
radii. The fact that particles of radii closer to r min are less 
probable affects the timescale for a successful adsorption 
to take place, since particles of larger radius cannot ad- 
sorb due to the overlap condition. The mean values of (77) 
are provided in Table |TJ Finally, the bulk case (a — > 00) 
has the highest skewness towards the smallest values of 
the radius; the mismatch between the distribution func- 
tions is larger (e.g., a = 1.7) at small radii for the inter- 
mediate values of the size dispersion. Geometrical con- 
straints imposed by small cell sizes effectively filters the 
smallest particles. The distribution function of the radii 
of adsorbed particles does show this trend with clearly 
asymmetric, non-Gaussian dependence on the normal- 
ized radius as show in Fig. [3j Moreover, at a cell size of 
a = 1.7, the distribution function of the adsorbed radii 
in the presence of a pattern closely follows that of regu- 
lar (non-patterned) substrate particularly for the highly 
polydisperse case. 

To further characterize the NICCA regime, we study 
the dependence of various quantities on the size of the 
cells, a. We first consider the jamming coverage, 9j, as 
a function of the cell size, a, for values of the size disper- 
sion ranging from 1% to 20% as shown in Fig. [4](a) . As 
observed in Sec. |III| values of the coverages, while within 
the SPCA regime, decay proportionally to (a + /3)~ 2 . 
Figure |4|a) shows the curve for the monodisperse case 
(offset down to 20% instead of 21.76% at a = 0.4 and 
(3 = 1.5 for clarity) as a guide for similar behavior in 
the polydisperse cases, and, of course, one observes that 
coverage values increase as the mean number of parti- 
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TABLE I. Values of the mean normalized radius, (77} for the 
various values of a = 1%, 5%, 10%, 15%, and 20% and the 
two values of a = 1.7 and the bulk case (a — > 00) belonging 
to the MPCA regime as used in Fig. [3] 





a = 1.7 


a — > 00 


1 


0.4851 


0.4792 


5 


0.4437 


0.4235 


10 


0.4007 


0.3686 


15 


0.3569 


0.3190 


20 


0.3109 


0.2712 



cles per cell is no longer unit. In the MPCA regime, 
an increase of the jamming coverage with the cell size is 
observed, as shown in Fig. [3]Ja) , regardless of the value 
of the size dispersion. For the highly monodisperse case 
of a size dispersion of 1%, the values of the jamming 
coverage with a closely follow those of the monodisperse 
case |17j , with features such as the cusp at the transition 
from up-to-two to up-to-three particles per cell present 
(a = (1 + \/3)/2 3 / 2 sa 0.96). More surprising is the 
presence of a faint signature of the cusp at a size disper- 
sion of 5% (see inset of Fig. [4f a)). Since the transition 
from SPCA to MPCA depends on the size of the small- 
est particle, an increase of the size dispersion leads to a 
lower transition threshold as shown in Fig. |4|a). This 
dependence can be observed in Fig. [2] for snapshots of 
the jammed state for size dispersions of 1% (Fig. |2^a)) 
and 20% (Fig.gb)) at a fixed a = (1 + \/3)/2 3 / 2 , which 
represents the transition from up-to-two to up-to-three 
particles per cell in a close packed arrangement in the 
monodisperse case. At this cell size and for a size dis- 
persion of 1%, only monomers and dimers are observed 
(Fig. [2ja)), though trimers could, in principle, be possi- 
ble. At a size dispersion of 20%, larger aggregates become 
possible, like pentamers (Fig. J2jb)). 

In general, the higher the values of the size dispersion, 
the higher the coverage, at the expense of a lower con- 
trol over both the number of particles adsorbed on each 
cell and their sizes. The mean particle number per cell 
is shown in Fig. [IJb). The mean adsorbed radius, 0, of 
the resulting film as a function of the cell size is shown 
in Fig. [4|c). The quantity reflects the balance between 
the population of small and large adsorbed particles as a 
function of the cell size and the degree of polydispersity of 
the particles. At small size dispersions, the values of the 
radii remain relatively undifferentiated with nearly sym- 
metric distributions as shown in Fig. [3](a) for the case of 
1%. As the size dispersion increases, the distribution of 
the adsorbed radii becomes more skewed towards smaller 
values. The overall result is reflected by a substantial de- 
crease of the average value of the mean adsorbed radius 
(at fixed a) as a function of the size dispersion as shown 
111 Fig. gc). This can be illustrated by considering the 



case of dimers for values of a above, though close, to 
a c . In this limit only small particles are allowed to ad- 
sorb in the formation of dimers, so the mean adsorbed 
radius must decrease. However, at a fixed value of a, 
say 0.8, and for increasing values of the size dispersion, 
the relative population of particles changes significantly 
as shown in Fig. [5] As the value of a c shifts to lower val- 
ues with increasing values of the size dispersion, dimers 
and trimers become possible. At the transitioning points, 
from monomer to dimer and from dimer to trimer, the 
cell acts as a filter to allow adsorption of the smallest par- 
ticles, i.e., particles with radii at or close to r min . As the 
values of a > a c increase the fraction of dimers grows, 
peaks, and decays; this feature reproduces itself on larger 
particle aggregates like trimers, tetramers, etc. The dis- 
tribution functions are skewed towards the large cell sizes. 
At large values of cell size (a > 1), although the mean 
number of particles per cell increases with the size disper- 
sion (Fig. |4](b)), yet the decrease in the mean adsorbed 
radius is such (Fig. (4^c)) that the coverage is no longer 
a monotonic increasing function of the size dispersion. 

Finally, we briefly address possible particle arrange- 
ments inside a cell. For example, dimers tend, on aver- 
age, to orient themselves along the diagonals of the cells. 
Considering cells sizes that accommodate up to three par- 
ticles, one can understand the underlying reasons. Par- 
ticles adsorbing at the center may block further particles 
to adsorb and become monomers. Now, if two particles 
adsorb along the diagonal, this blocks the possibility of 
a third particle to adsorb, thus becoming a dimer. Fi- 
nally, two particles adsorbed along an edge permit a third 
particle to fit in and form a trimer. 



V. CONCLUDING REMARKS 

We performed extensive Monte Carlo simulations of 
the irreversible adsorption of polydisperse disks on a pat- 
terned substrate. A pattern consisting of square cells 
positioned at the vertices of a square lattice was consid- 
ered, but extensions and generalizations to other lattice 
arrangements and cell shapes are straightforward. We 
used a physically motivated truncated Gaussian-size dis- 
tribution to model the polydispersity of the particles at- 
tempting adsorption and size dispersions up-to 20% of 
the average radius. The model extends a previous study 
of the irreversible adsorption of monodispersed particles 
on patterned substrates [17] ■ In the present work, we 
focused on the noninteracting cell-cell adsorption regime 
(NICCA) and on the jammed state properties. 

The model is suitable to describe relevant features of 
colloidal particle adsorption under the assumption that 
no particle can fit under another particle of larger ra- 
dius as discussed in Sec. [TTJ Even though the modeling 
relies on the excluded volume interaction, several quan- 
tities are generally valid, even if more realistic interac- 
tions are taken into account, like Coulomb and van der 
Waals. For example, the various transition points of the 
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a 

FIG. 5. (Color online) Density of monomers (red, solid 
line), dimers (green, dashed line), trimers (blue, dotted line), 
tetramers (magenta, fine dot line), and pentamers (light blue, 
dash-dot line) for different values of the size dispersion: a) 1%, 
b) 5%, c) 10%, d) 15%, and e) 20%. 



phase diagram, as, e.g., the transition from monomer to 
dimer solely depends on geometrical parameters such as 
the minimum and maximum values of the radii of parti- 
cles attempting adsorption. Hence, these transition val- 
ues are not affected by the inclusion of more complex 
interactions especially for the more interesting case of 
small cell sizes, where geometrical constraints are more 
significant and one observes substantial departure from 
the bulk (regular) substrates. In contrast, the interparti- 
cle distance of the adsorbed particles will depend on the 
choice of interaction potential, or the presence of hydro- 
dynamic effects, since the quantity depends on the way 
adsorbed particles interact with a landing particle. Cov- 
erage represents another quantity that depends on the 
interactions taken into account. 

Coverage efficiency, for the excluded volume interac- 
tion case, improves with size dispersion even accounting 
for the fact that the adsorption of small particles is fa- 
vored. The presence of a pattern favors adsorption of 
small particles for cell sizes near, but above, the transi- 
tion points, e.g., from monomer to dimer and dimer to 
trimer. 

Finally, the present study provides insight on the val- 
ues of the parameters required to tune the average num- 
ber of particles adsorbed per cell. 



ACKNOWLEDGMENTS 



A.C. thanks both Fundacao para a Ciencia e a Tecnolo- 
gia (SFRH/BPD/34375/2007) and Fundacao Calouste 
Gulbenkian for fellowships to visit Los Alamos National 
Laboratory and acknowledges the warm hospitality of 
the T-l Group at Los Alamos National Laboratory. We 
thank Vladimir Privman and Cynthia Reichhardt for use- 
ful comments on the manuscript. 



[1] E. Kumacheva, R. K. Golding, M. Allard, and E. H. 
Sargent, Adv. Mater., 14, 221 (2002). 

[2] C.-A. Fustin, G. Glasser, H. W. Spiess, and U. Jonas, 
Adv. Mater., 15, 1025 (2003). 

[3] G. Zhu, C. Wang, Y. Zhang, N. Guo, Y. Zhao, R. Wang, 
S. Qiu, Y. Wei, and R. H. Baughman, Chem. Eur. J., 
10, 4750 (2004). 

[4] C. Burda, X. Chen, R. Narayanan, and M. A. El-Sayed, 
Chem. Rev., 105, 1025 (2005). 

[5] P. C. Lewis, E. Kumacheva, M. Allard, and E. H. Sar- 
gent, J. Disp. Sci. and Technol., 26, 259 (2005). 



[6] E. O'Connor, A. ORiordan, H. Doyle, S. Moynihan, 
A. Cuddihy, and G. Redmond, App. Phys. Lett., 86, 
201114 (2005). 

[7] L. Yang, Z. Yang, W. Cao, L. Chen, J. Xu, and H. Zhang, 

J. Phys. Chem. B, 109, 11501 (2005). 
[8] N. L. Jeon, W. Lin, M. K. Erhardt, G. S. Girolami, and 

R. G. Nuzzo, Langmuir, 13, 3833 (1997). 
[9] J. Y. Chen, J. F. Klemic, and M. Elimelech, Nano Lett., 

2, 393 (2002). 

[10] M. Elimelech, J. Y. Chen, and Z. A. Kuznar, Langmuir, 
19, 6594 (2003). 



9 



[11] J. Joo, B. Y. Chow, and J. M. Jacobson, Nano Lett., 6, 
2021 (2006). 

[12] D. Xia, D. Li, Y. Luo, and S. R. J. Brueck, Adv. Mat., 
18, 930 (2006). 

[13] K. M. Chen, X. Jiang, L. C. Kimerling, and P. T. Ham- 
mond, Langmuir, 16, 7825 (2000). 

[14] H. Zheng, M. F. Rubner, and P. T. Hammond, Lang- 
muir, 18, 4505 (2002). 

[15] A. Cadilhe, N. A. M. Araiijo, and V. Privman, J. Phys.: 
Condens. Matter, 19, 065124 (2007). 

[16] Z. Adamczyk, J. Barbasz, and M. Nattich, Langmuir, 
24, 1756 (2008). 

[17] N. A. M. Araiijo, A. Cadilhe, and V. Privman, Phys. 
Rev. E, 77, 031603 (2008). 

[18] P. Meakin and R. Jullien, Phys. Rev. A, 46, 2029 (1992). 

[19] P. Meakin and R. Jullien, Physica A, 187, 475 (1992). 

[20] Z. Adamczyk, B. Siwek, M. Zembala, and P. Werohski, 
J. Coll. Interf. Sci., 185, 236 (1997). 

[21] O. Gromenko and V. Privman, Phys. Rev. E, 79, 011104 
(2009). 

[22] O. Gromenko and V. Privman, Phys. Rev. E, 79, 042103 

(2009) . 

[23] J. Bafaluy, B. Senger, J.-C. Voegel, and P. Schaaf, Phys. 
Rev. Lett., 70, 623 (1993). 

[24] B. Senger, J.-C. Voegel, P. Schaaf, A. Johner, A. Schmitt, 
and J. Talbot, Phys. Rev. A, 44, 6926 (1991). 

[25] B. Senger, P. Schaaf, J. C. Voegel, A. Johner, A. Schmitt, 
and J. Talbot, J. Chem. Phys., 97, 3813 (1992). 

[26] B. Senger, J. Talbot, P. Schaaf, A. Schmitt, and J.-C. 
Voegel, Europhys. Lett., 21, 135 (1993). 

[27] E. Kumacheva, P. Garstecki, H. Wu, and G. M. White- 
sides, Phys. Rev. Lett., 91, 128301 (2003). 

[28] R. Duffadar, S. Kalasin, J. M. Davis, and M. M. Santore, 
J. Coll. Interf. Sci., 337, 396 (2009). 

[29] V. P. Zhdanov and B. Kasemo, Coll. Surf. B, 76, 28 

(2010) . 

[30] This might not be necessarily true for certain systems 
involving proteins as discussed, e.g., [33U57] . 

[31] J. W. Evans, Rev. Mod. Phys., 65, 1281 (1993). 

[32] V. Privman, J. Adhesion, 74, 421 (2000). 

[33] V. Privman, Coll. and Surf. A, 165, 231 (2000). 

[34] V. Privman (ed.), Adhesion of submicron particles on 
solid surfaces, Coll. and Surfaces A, 165, Nos. 1-3 
(2000). 

[35] J. Talbot, G. Tarjus, P. R. van Tassel, and P. Viot, Coll. 



and Surf. A, 165, 287 (2000). 
[36] P. Schaaf, J.-C. Voegel, and B. Senger, J. Phys. Chem 

B, 104, 2204 (2000). 
[37] T. W. Odom, V. R. Thalladi, J. C. Love, and G. M. 

Whitesides, J. Am. Chem. Soc, 124, 12112 (2002). 
[38] K. E. Paul, M. Prentiss, and G. M. Whitesides, Adv. 

Func. Mater., 13, 259 (2003). 
[39] S. M. de Sousa Pereira, M. A. Martins, T. Trindade, I. M. 

Watson, D. Zhu, and C. J. Humphreys, Adv. Mater., 20, 

1038 (2008). 

[40] Z. Adamczyk and P. Werohski, J. Chem. Phys., 108, 9851 
(1998). 

[41] A. Cadilhe and V. Privman, Mod. Phys. Lett. B, 18, 207 
(2004). 

[42] J. C. Voegel, N. de Baillou, and A. Schmitt, Coll. Surf., 
16, 289 (1985). 

[43] E. Pefferkorn, A. Carroy, and R. Varoqui, J. Polymer 
Sci. Polymer. Phys. Ed., 23, 1997 (1985). 

[44] H. Nygren, S. Alaeddin, I. Lundstrom, and K.-E. Mag- 
nusson, Biophys. Chem., 49, 263 (1994). 

[45] T. F. Tadros and J. Gregory (ed.), Particle Deposition at 
the Solid-Liquid Interface, Colloids Surf., 39 (1989). 

[46] M. Elimelech, J. Gregory, X. Jia, and R. A. Williams, 
Particle Deposition & Aggregation. Measurement, Mod- 
eling and Simulation (Butterworth-Heinemann Woburn, 
MA, 1995). 

[47] D. S. Rimai and L. H. Sharpe, eds., Advances in Particle 
Adhesion (Gordon and Breach Publishers, Amsterdam, 
1996). 

[48] J. Feder, J. Theor. Biology, 87, 237 (1980). 
[49] J. Talbot and P. Schaaf, Phys. Rev. A, 40, 422 (1989). 
[50] The fact we take ±2<r represents just a choice, not a 
requirement. 

[51] This is a direct consequence of considering symmetric 

values of a around the mean value. 
[52] The actual definition of 7 is provided in the next section. 
[53] Further discussions can be found in [18 20 . 
[54] N. A. M. Araiijo and A. Cadilhe, Unpublished. 
[55] V. Privman, J.-S. Wang, and P. Nielaba, Phys. Rev. B, 

43, 3366 (1991). 
[56] J.-S. Wang, Int. J. Mod. Phys. C, 5, 707 (1994). 
[57] V. Privman, H. L. Frisch, N. Ryde, and E. Matijevic, J. 

Chem. Soc. Faraday Trans., 87, 1371 (1991). 



